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Abstract. Pulsars are powerful sources of radiation across the electromagnetic spectrum. This 
paper highlights some theoretical insights into non-thermal, magnetospheric pulsar gamma-ray 
radiation. These advances have been driven by NASA's Fermi mission, launched in mid-2008. The 
Large Area Telescope (LAT) instrument on Fermi has afforded the discrimination between polar cap 
and slot gap/outer gap acceleration zones in young and middle-aged pulsars. Altitude discernment 
using the highest energy pulsar photons will be addressed, as will spectroscopic interpretation of 
the primary radiation mechanism in the LAT band, connecting to both polar cap/slot gap and outer 
gap scenarios. Focuses will mostly be on curvature radiation and magnetic pair creation, including 
population trends that may afford probes of the magnetospheric accelerating potential. 
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INTRODUCTION 

Prior to the launch of the Compton Gamma-Ray Observatory (CGRO) in 1991, the 
gamma-ray pulsar database was extremely limited. CGRO extended the population to 
seven sources, all of them young pulsars, and only Geminga had no definitive observa- 
tion of a radio pulsar counterpart. There were also a few marginal CGRO detections, 
including the millisecond pulsar PSR J02 18-1-4232. Following the launch of the Fermi 
Gamma- Ray Space Telescope on June 11, 2008, the gamma-ray pulsar database has 
grown rapidly, with over 75 sources now observed with high significance. Most of these 
are again young pulsars with higher spin-down luminosities than the average for ra- 
dio pulsars. Yet, a sizable millisecond pulsar (MSP) population has emerged, at present 
numbering more than twenty. Moreover, the improved photon counts captured by the 
Fermi Large Area Telescope has enabled blind pulsation searches, those without prior 
knowledge of a radio timing ephemeris. The LAT has discovered over 20 pulsars in blind 
searches of the gamma-ray data alone, and only a few have since had radio confirmation. 

The improved pulse-phase spectroscopy enabled by Fermi has clarified the probable 
locale of the gamma-ray emission region. The two major competing models for the site 
of this emission are the polar cap (PC) near the stellar surface (e.g. [10, 12]) and the 
outer gap (OG) near the light cylinder (e.g. [9, 20]). Super-exponential turnovers due to 
magnetic pair creation are expected for PC models at the maximum photon energies 
in the GeV band, as opposed to exponential shapes for OG models. Discrimination 
by Fermi-LP^I between these two shapes was anticipated prior to launch [19]. In the 
vast majority of LAT pulsar spectra, the turnovers are exponential, thereby favoring the 
OG scenario [2]. Furthermore, the separation of peaks in pulsars with more than one 
pulse peak is generally of the order of 0.3-0.5 in phase, indicating a broad radiation 



cone that is consistent only with high ahitude emission locales for dipolar fields (e.g. 
[26]). This evidence has ushered in a new era for the y-ray pulsar paradigm. Here, 
some perspectives pertaining to this new era are offered, focusing first on lower bounds 
to the emission radius from magnetic pair creation physics. In addition, the GeV-band 
turnover energies ^max are assessed in terms of the canonical radiation-reaction-limited 
curvature emission picture, touching upon population characteristics for y-ray pulsars. 

MAGNETIC PAIR CREATION BOUNDS FOR THE MINIMUM 

ALTITUDE OF EMISSION 

The emergence of gamma-rays from pulsars is probably intimately coupled to the gen- 
eration of electron-positron pairs in their magnetospheres. Such a connection between 
radiation emission and pair creation was posited in early radio-pulsar pictures [21, 24]. 
In the inner magnetosphere of young pulsars, the dominant pair production process in- 
volves a single photon interacting with the B field, 7 — e+e^ , with its non-magnetic, 
quantum counterpart 77— ?■ e^e^ only becoming competitive at high altitudes. Hence, 
7—7' e+e^ is the preserve of polar cap (PC; [10]) and slot gap (SG; [18]) gamma-ray 
pulsar models. It is a first-order QED process that is kinematically forbidden in field-free 
regions, but can take place in an external magnetic field, which can absorb momentum 
perpendicular to B. Its generic physics properties are summarized in [6]. Due to energy 
conservation, its rate 7?pp possesses an absolute threshold energy, Cy > 2/ sin 0)^3 . Here 
01jB is the angle the photon momentum vector k makes with the magnetic field B, and 
the photon energy is in units of m^c^ , a convention adopted hereafter. 

Above threshold, EP^ , which is averaged here over photon polarizations, exhibits a 
large number of resonances, corresponding to thresholds for the creation of pairs in 
excited states; these aggregate to exhibit a characteristic sawtooth structure (e.g. [6, 11]). 
In pulsar contexts, considerable ranges of field strengths and photon angles 9^ are 
sampled during magneto spheric propagation: the resulting convolution smears out the 
sawtooth appearance into a continuum. In this regime, the dependence on £y can be 
approximated by a compact asymptotic expression (e.g. [13]; and references in [6]): 

^■"'^ ~ ly^^^ sin0kB exp|-^| , T = -^CysinekB < 1 , (1) 

where T is the critical asymptotic expansion parameter. Here, Uf = j (fic) is the 
fine structure constant. Arc = h/{mec) = 3.862 x 10^^^ cm is the electron Compton 
wavelength over In, and 5cr = m^c^ / (eh) = 4.413 x 10^^ Gauss is the quantum critical 
field, where the cyclotron energy equals mgC^ . Near threshold, this asymptotic result 
is an imprecise description of 7?pp [11], and can be improved considerably by careful 
treatment of the associated mildly-relativistic regimes for the produced pairs (e.g. [4]). 

In polar cap pulsar models [10, 21, 24], primary curvature photon emission is emitted 
at very small angles to the magnetic field, well below pair threshold (e.g. [8]). These 
photons will convert into pairs only after traveling a distance s that is a fraction of the 
field line radius of curvature pc . Above the polar cap, the radius of field curvature is 
pc = [Prc/2nY/^ for a pulsar period P, and exceeds the neutron star radius i?Ns . Pair 



creation ensues for small angles such that sin^j^B ^ s/pc, and the argument T of the 
exponential becomes a fraction of unity, i.e., when £y5sin0kB ^ 0.25cr. Accordingly, 
in young pulsars, pair production will usually occur somewhat or well above threshold 
when B <^ 0.15cr, and the attenuation mean free path App ~ pc/ £ max{2, Q.2Bcx/B} 
will then be much less than R^^ . Consequently, the optical depth Ty(£y) for one-photon 
pair opacity will be a rapidly-varying function of Cy, producing a dramatic attenuation 
signature at the escape energy e^^^ , which typically lies above 100 MeV (see [8]). 

The upshot is that attenuation by 7— e+e^ should impose a super-exponential cutoff 
in pulsar 7-ray spectra of the approximate form dny/d£y°^ exp{ — aexp[— Ccsc/Cy]} for 
some constant a . If Cosc ~ £max = £^MAx/"^eC^ , the severity of the cutoff should produce 
[10, 12] a sharp turnover that is distinguishable from those generated by other mecha- 
nisms. It was widely anticipated [19] that the Fermi-LP^ instrument would discern the 
precise shape of such cutoffs in Vela and other pulsars, and accordingly discriminate 
between the polar cap and outer gap models for their high energy emission. It has in fact 
done so. Early on in the first year of observations by Fermi, the high count rates for Vela 
permitted the exclusion of super-exponential turnovers at a high level of significance 
[1]. This has become possible for the pool of 39 young pulsars listed in the Ferm/-LAT 
Pulsar Catalog in [2], all of which reveal significant turnovers in the 1-10 GeV band 
consistent with simple exponentials. This clearly demonstrates that 7 — e+e^ is not 
producing the observed turnovers, proving that Ty(ey) < 1 for Ey < Emax ■ 

This undeniable absence of the signature of active magnetic pair creation in gamma- 
ray pulsar signals can be inverted to provide robust lower bounds to emission altitudes, 
a protocol that has been adopted in a number of Fermi-LAI pulsar papers. Defining a 
marginal magnetic pair opacity criterion Ty(eMAx) ~ 1 near the pulsar surface, rearrange- 
ment yields an approximate dependence of pair creation cutoff energies Emax on Bq , Rq 
and pulsar period P (in seconds). Setting ^max = ^MAx'^ef ^ , this can be summarized in 
a relation [5] that corresponds to near-surface emission (r < 2i?Ns ): 

1 /2 f 3 ^ 

£max ~ OAVp maxi 1, i GeV . (2) 

This encapsulates the character of accurate numerics derived from codes employed in 
[14] and [8], which include pronounced effects of general relativity on spacetime curva- 
ture, field enhancement and photon energy. In flat spacetime, the altitudinal dependence 
is weaker, corresponding to Emax °^ (r/R^s)^^'^ (e-g- [16, 27]) above the magnetic poles. 

The pair opacity Cmax trend can be expressed alternatively as a combination of 
observables appropriate for this inner gap (polar cap) discussion. The form is 

logio^l^ ^ a:c + ^logioZiG+0.25 , = -Iogio5i2 + 5 logio/' • (3) 

which clearly encapsulates an anti-correlation with P, as does Eq. (2), since Bq = 
3.2 X lO^'^{PPy^-^ [17]. Here Xig = r/Rt^s is the scaling of the radius appropriate for 
polar cap considerations. Figure 1 plots the phase space corresponding to Eq. (3), using 
the measured spectroscopic cutoffs Ec for young Ferm/-LAT pulsars in Table 4 of [2] to 



Magnetic Pair Creation Cutoffs in Pulsars 
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FIGURE 1. The inferred maximum pulsar emission energies ^max as a function of the key polar cap 
(inner gap) magnetic pair creation attenuation parameter Giq defined in Eq. (3). The spectral parameters 
used are the cutoff energies E^. listed in Table 4 of [2], the First Ferm;-LAT Catalog of Gamma-ray 
Pulsars. For the 37 pulsars constituting the black points, Emax = 2.5Ec was set; for the Crab and Vela, 
Emax = 3-5Ec (see text). The high-field CGRO Comptel pulsar B 1509-58 (J15 13-5908) is also displayed. 
No LAT millisecond pulsars are depicted. The fiducial linear relationship in Eq. (3) is displayed as the 
diagonal lines for four different emission altitudes, r/Rj^s = 1 (surface), 2,3,4, as labelled. 



benchmark the maximum energies ^max • Since these cutoffs are exponential in charac- 
ter, the action of magnetic pair creation can be ruled out at energies below the window 
around Ec. For this reason, a representative setting of ^max = 2.5£'c was deployed in 
the Figure for most of the pulsars, the set of 37 with points in black. This value can be 
adjusted on a case-by-case basis, depending on the photon counting statistics. Accord- 
ingly, the alternative value of ^max = S-SE^- was adopted for the Crab and Vela pulsars. 
Motivations for these exceptions are as follows: the Crab, possesses high energy emis- 
sion reflected in its pulsed detection by MAGIC [3], and Vela, because its impressive 
count rate as the brightest LAT-band pulsar in the sky permits extraordinary spectro- 
scopic determination of its turnover [1]. The uncertainties in E^- were scaled similarly, 
yielding the ^max error bars in the Figure. PSR B 1509-58, the young, high-field pulsar 
detected by Comptel on CGRO, but not by EGRET, was also placed in the Figure using 
an inferred ^max value from the spectral results displayed in [14]. It is the only gamma- 
ray pulsar that might be consistent with 7 — )■ e+e^ attenuation operating near the stellar 
surface. Observe that no millisecond pulsars (with Oiq < —2) are exhibited in Fig. 1 as 
their surface fields are so low that 7 — pair creation is exponentially improbable in 
their magnetospheres. The phase space plot in Figure 3 clearly demonstrates that there 
is no correlation between £max 

and OiG or P in the LAT pulsar population. 



CURVATURE RADIATION AT HIGH ALTITUDES 



Given the above discussion it is evident that the Fermi GeV-band spectral cutoffs in 
young pulsars are probably related to the primary emission mechanism. Dating from the 
early models of [24], this mechanism is postulated to be curvature radiation, the ele- 
mentary dipole radiation resulting from pairs accelerating in the curved magnetospheric 
fields. The formalism for the curvature radiation emissivity is identical to that for classi- 
cal synchrotron radiation [15, 22], but with the electron Larmor radius r^ = YetrieC^/ (eB) 
being replaced by the radius of field curvature, pc . Accordingly, the characteristic di- 
mensionless (i.e., in units of mgC^ ) photon energy scales for curvature radiation, £c , and 
synchrotron emission , are 

3^A- 3^Ar 35? 

^^■^2^p; ' = 2^7f ^^"^^ = 2 5^^^^"^" • ^"^^ 

For curvature emission to seed pulsar gamma-ray signals, £c must exceed a few GeV in 
the magnetosphere; this is routinely satisfied in polar cap (e.g. [12]), slot gap [18] and 
outer gap (e.g. [20]) models, for both young pulsars and millisecond ones. 

At its characteristic photon energy, the differential emissivity per electron scales as 
afc/^cYi^ (^-S- Eq. (2) of [7]) for both the curvature and synchrotron processes. 
When this emissivity scale is multiplied by the square of the photon energy scale 
in Eq. (4), the corresponding magnitudes for the curvature radiation and synchrotron 
cooling rates result; in full form these rates are 

Pc fg 

Here rg = meC^/{eB) — ^c^cr/B is c divided by the cyclotron frequency, and ro = 
OCfXc = e^ I {me(?) is the classical electron radius. The precise numerical factors can be 
derived using the synchrotron formulation of [22]; for the curvature process this requires 
the specialization to pitch angles Qe = n/l. From Eq. (5) one derives the criterion for 
the dominance of a synchrotron signal over a curvature one: Je < pc sin 0e/rg . While 
0^ ^ is presumed for the primary emission, pair creation usually generates secondary 
photons at significant angles to the field, so that this inequality is almost always satisfied: 
synchrotron photons therefore abound in subsequent generations of a cascade. Since 
Pc/i'g ^ 10^^ near the polar cap in young pulsars, only very small Og are required to 
spawn the swamping of primary curvature emission by a synchrotron component in the 
inner magnetosphere. This possibility should be borne in mind for future refinements of 
the gamma-ray pulsar paradigm. 

The maximum Lorentz factor of the electrons emitting in gamma-ray pulsars is a 
function of the accelerating potential adopted in a given model. If the electric field com- 
ponent along the local magnetic field is (presumed positive), then the energy gain 

rate in the electrostatic gap is %ccineC^ = eE\^\c, i.e. the electron speed c times the force. 
The maximum possible electron Lorentz factor /max is realized when this gain equals 
the loss rate due to radiative cooling, which is the radiation-reaction limit generally as- 
cribed to curvature emission in gamma-ray pulsar models. Setting —eE\\ / (mgc) equal to 



the first result in Eq. (5) yields a maximum Lorentz factor /max ~ [3Pc£^||/(2e)] , 
which can then be inserted into Eq. (4) to specify the approximate maximum photon 
energy (in units of mgC^ ) for curvature radiation reaction-limited acceleration: 

This form for Cmax is routinely cited in Fermz-LAT publications on gamma-ray pulsars 
(e.g. [2]). In principal, the observed cutoff energy can be less than this value if atten- 
uation by pair production, magnetic or two-photon, is efficient, or the acceleration is 
terminated by mechanisms other than curvature cooling. Hence, the observed GeV-band 
turnovers in dozens of young pulsars in the First Ferm?-LAT Pulsar Catalog provide 
lower bounds to the accelerating potential, modulo the altitude of emission. 

A natural, fiducial scale for the parallel electric field can be obtained from Lorentz 
transformations associated with a rotating magneto sphere. This establishes rQ.^B/c 
as the co-rotational component of E , providing the ideal MHD contribution. Yet within 
the plasmasphere, a non-co-rotational, parallel component Fy emerges due to departures 
from Goldreich- Julian current flow (e.g. [23, 25]): it is putatively of the same order of 
scaling ~ rQ.B/c. This is the most optimistic scenario, being adopted for example in 
the seminal outer gap pulsar model of [9], and leads to E\\ ~ e/ (ro-^c)Z^^Lc/^cr for 
X — r/Ri^c describing the scaling of the acceleration/emission radius in terms of the 
light cylinder radius 7?lc = Pc/ (In) that will be employed hereafter. These define so- 
called thick outer gaps (e.g. [28]). Here 5lc is the field strength at r = i?Lc > essentially 
along the last open field line. Note that the dependence of Fy , B^^ and the radius of 
curvature on the obliquity of the rotator, and colatitude 0, is omitted in this discussion. 
Using the scaled radius of curvature Vc = Pc/ (z^lc) » one can then insert this nominal 
value for Fy into Eq. (6). This can be cast in terms of observables P and P. For a 

surface polar field Bq = 3.2 x 10^^ (PP) = 10^^5i2 Gauss and a neutron star radius of 
i?Ns = 10^ cm, since = (Rns/Rlc)^ , setting P = OAP i sec yields a result of 

£max ^ ^^(£||5i2)3/^(/'-i)-^/4GeV . (7) 

This scaling is commensurate with that obtained in Eq. (4) of [20]. Clearly Fmax 
rises monotonically as the acceleration locale moves from the light cylinder towards 
the neutron star surface, due to the rapidly rising inductive 5 -field. An additional 
electrostatic decrement factor e | < 1 has been introduced in Eq. (7) to facilitate the 
subsequent discussion. It is an acceleration efficiency factor, representing the departure 
of £■[! below the fiducial, optimistic value of E\\ ~ rQ.B/c ; in general it is a function of 
the altitude parameter x > the colatitude and the rotator obliquity angle a . 

To interpret this outer magneto spheric curvature radiation cutoff scaling. Figure 2 dis- 
plays the phase space corresponding to Eq. (7), plotting logjo^c versus the appropriate 
combination of observables P and P; in this graph, the cutoff energy E^- serves as a 
proxy for Fmax • For the curvature radiation expectation in an outer gap (OG) near the 



Curvature Radiation Cutoffs in Fermi -LAT Pulsars 
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FIGURE 2. Observed Fermi-LAJ^ pulsar emission cutoff energies as a function of the key outer gap 
curvature radiation reaction-limited acceleration parameter Gqc defined in Eq. (8). The data come 
from the First Ferm;-LAT Pulsar Catalog [2], with classic spinning-down pulsars (SDPs) in black, and 
millisecond pulsars (MSPs) in blue, as indicated. The Crab pulsar point from the Fermi dataset is 
highlighted, noting that its Ec from the MAGIC detection [3] is considerably higher The diagonal lines 
labelled with x = '"/^lc values constitute solutions of Eq. (8) for en = 1 and = 1 . 



light cylinder, this is (for jj^ = ^^'^ ) 

E 3 7 

^"^lorSiV ^°°~^°^io^"^-^^ ' = 4logio5i2-4logioi' . (8) 

The graph is derived from the LAT pulsar data in Tables 1 and 4 of [2]. Since the 
Emax values appear in a relatively narrow range, with considerable spread in Oqg, 
there appears to be no clear linear correlation, just the suggestion of one. The diagonal 
lines corresponding to x = 1,0.1,10^^ in the Figure provide an approximate guide 
to read off the scaled altitude. It is clear that this case of maximal cannot quite 
accommodate the turnover energies of the LAT pulsars being attributed to cooling- 
limited curvature emission within the light cylinder: the observations mandate that the 
electrostatic decrement lie in the range 0.05 ^ £|| < 1 if X = ^ ■ This essentially implies 
acceleration potentials somewhat weaker than the co-rotational benchmark rQ.B/c, or 
considerably so (i.e., e|| < 10^^) if the outer gap emission altitude lies well within 
the light cylinder radius, i.e. x ^0.1. This observational diagnostic becomes less 
constraining if some two-photon pair attenuation is active in the emission region, thereby 
lowering the observed cutoff energy Ec below the curvature radiation ^max • 



The phase space plot in Fig. 2 exhibits a remarkable property — that the youngish, 
dipole-torque spin-down pulsars (SDPs) and old, recycled millisecond pulsars (MSPs) 
occupy the same locale. The .)c-coordinate, cTog ^ is an approximate linear function of 
logio^Lc , and the field 5lc is fairly similar for the two populations (see Table 1 of [2], 
and compare Fig. 2 here with their Figure 7 that displays the array of Fermz-LAT pulsar 
£max values versus logjo^Lc )• This approximate coincidence of 5lc values in SDPs and 
MSPs was, of course, well-known for the radio pulsar population. Yet, the clustering of 
the gamma-ray Emax around similar values for SDPs and MSPs is a new insight enabled 
by Fermi. It is strongly suggestive that similar physical acceleration locales and radiative 
emission mechanisms operate in the two populations. 

Concluding Remarks: The picture of high altitudes for y-ray emission in young pulsars 
is now well-established. Yet questions still remain. Is the primary mechanism curvature 
emission, or can small-angle synchrotron radiation contribute significantly to the signal? 
What is the role of two-photon pair creation in controlling the range of observed Cmax ? 
Does magnetic pair production at low altitudes still enhance the multiplicity of pairs 
being pumped into surrounding pulsar wind nebulae? These issues will capture the focus 
of theorists in the years to come as the Fermi pulsar legacy further unfolds. 
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